We present a mode identification based on new high- 
resolution time-series spectra of the non-radially pulsat- 
ing 5 Scuti star FG Vir (HD 106384, V = 6.57, A5V). 
From 2002 February to June a global Delta Scuti Network 
(DSN) campaign, utilizing high-resolution spectroscopy 
and simultaneous photometry has been conducted for 
FG Vir in order to provide a theoretical pulsation model. 
In this campaign we have acquired 969 Echelle spectra 
covering 147 hours at six observatories. The mode iden- 
tification was carried out by analyzing line profile varia- 
' tions by means of the Fourier parameter fit method, where 
f^i , the observational Fourier parameters across the line are 
' fitted with theoretical values. This method is especially 
well suited for determining the azimuthal order m of non- 
£^ radial pulsation modes and thus complementary with the 
^ | method of Daszynska-Daszkiewicz (2002) which does best 
at identifying the degree £.15 frequencies between 9.2 and 
33.5 d _1 were detected spectroscopically. We determined 
the azimuthal order m of 12 modes and constrained their 
harmonic degree I. Only modes of low degree (I < 4) 

■ were detected, most of them having axisymmetric charac- 
l/"") \ ter mainly due to the relatively low projected rotational 

velocity of FG Vir. The detected non-axisymmetric modes 
have azimuthal orders between -2 and 1. We derived an 

■ inclination of 19°, which implies an equatorial rotational 
\Q ' rate of 66 kms . 

o \ 
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Abstract. 



1. Introduction 

Asteroseismology permits the probing of the theory of the 
interior structure and evolution of oscillating stars from 
the precise measurement of their light variations. 5 Scuti 
stars are particularly useful for this technique due to their 
large number of simultaneously excited pulsation modes 
having periods of a few hours. For these stars, every single 
detected pulsation frequency provides new independent 
physical information, allowing us to distinguish between 
different theoretical models. 

The techniques of asteroseismology rely on a detailed 
comparison between the observed and modeled pulsation 
frequencies. In order to reduce the number of free parame- 
ters, we need to identify the pulsational quantum numbers 
associated with each mode for as many frequencies as pos- 
sible. High resolution time-series spectroscopy of absorp- 
tion lines is a very effective tool for mode identification 
and permits not only the typing of I and m, but can also 
set constraints to the stellar inclination, the intrinsic pul- 
sation amplitude and involved temperature variations. 

A key object for successful asteroseismic studies is FG 
Vir, one of the best studied S Scuti stars. It has been the 
target of numerous photometric and spectroscopic observ- 
ing campaigns during the past 20 years. Despite intense ef- 
forts no definitive pulsational modeling has been achieved 
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yet, although we already have detailed knowledge about 
its rich frequency spectrum from the observational side. 
This is mainly due to a lack of good identifications of its 
pulsation modes. The determination of £ and to for at 
least a few modes would decrease the number of degrees 
of freedom of possible theoretical pulsation models. 

We applied the Fourier parameter fit (FPF) method 
(Zima 2006), hereafter referred to as Paper I, to spectra 
of FG Vir for the identification of the detected pulsation 
modes. This method is ideally suited for narrow lined stars 
pulsating with low-degree modes and enables the determi- 
nation of £ and to, as well as other pulsational parameters. 

This paper is constructed in the following way. First, 
we provide a short summary of the observational history of 
FG Vir, followed by a description of the new spectroscopic 
measurements and the data reduction. The spectroscopic 
detection of pulsation frequencies from the analysis of line 
profiles is discussed. Finally, we report the mode identifi- 
cation by means of the FPF method. 



2. Previous spectroscopic observations of FG Vir 

Eggcn (1971) discovered the variability of FG Vir in 
1970 during three hours of photometry detecting one pe- 
riod of 0.07 d (14.2 d _1 ) and a semi-amplitude of 25 
mmag. Subsequent photometric observations of this ob- 
ject (Dawson et al. 1995) revealed a rich pulsation spec- 
trum making it a target of further extended studies in the 
following years. 
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Pioneering work has been done by Mantegazza et al. 
(1994), who realized that only a combined photomet- 
ric and spectroscopic approach would yield significant 
progress for pulsational modeling. Eight nights of single- 
site differential photometry in V revealed seven frequen- 
cies. The high-resolution spectra obtained during a single 
night showed that FG Vir has a low projected rotational 
velocity of about vsini=21 kms -1 resulting in narrow ab- 
sorption lines. The dominant mode was identified as prob- 
ably the first radial overtone mode - contrary to the results 
of Dawson et al. (1995). 

Viskum et al. (1998) applied the equivalent width 
method to FG Vir in order to derive t for eight frequencies. 
Their results were in good agreement with the identifica- 
tion by Breger et al. (1999), who applied a photometric 
mode identification based on phase shifts and amplitude 
ratios of frequencies in two different pass bands. Balona 
(2000) showed that the two applied techniques are not in- 
dependent. Furthermore, Viskum et al. (1998) claimed a 
rotation period of 3.5 d derived from variations in the time 
series which would indicate an equatorial rotation velocity 
of about 33 kms -1 and an inclination of 40°. 

A detailed atmospheric study of FG Vir was made by 
Mittermayer & Weiss (2003) in order to derive its chemical 
composition. Their abundance analysis showed a solar-like 
pattern with slight overabundances of Ba and an under- 
abundancc of C and S. A metal overabundance claimed by 
Russell (1995) could not be confirmed. The analysis was 
hampered by the pulsational broadening of the absorp- 
tion lines, which required an unrealistically high macro- 
turbulence in order to model the spectra. 

Recently, Mantegazza & Poretti (2002) published re- 
sults from five nights of high-resolution time-series spec- 
troscopy obtained during the 1995 DSN campaign of FG 
Vir. They analyzed the line profile variations with the 
pixel-by-pixcl method and were able to detect ten pre- 
viously known frequencies in the radial velocity and pixel- 
by-pixel variations. Their identification revealed eight ax- 
isymmetric modes (m = 0) and only two modes with 
i > 1. They derived an inclination value of 15° which 
is much lower than the value proposed by Viskum et al. 
(1998), leading to a much higher equatorial rotation ve- 
locity of 80 kms -1 . By combining flux and velocity mea- 
surements they were able to derive the following physi- 
cal parameters for the dominant mode: £ = 1, m = 0, 
a s = 0.00147 ± 0.026, i = 15°, f r = -4.3 and /, = 14, 
where f r and /j are the real and imaginary part of /, the 
ratio of the stellar flux to radius variations. 

Meanwhile, a number of extended global photomet- 
ric observing campaigns, carried out by the Delta Scuti 
Network (Breger et al. 1998, 2004, 2005) revealed 79 pul- 
sation and combination frequencies. The existence of sev- 
eral very close frequencies, some being separated less than 
0.01 d _1 (Breger & Pamyatnykh 2006), makes seismic 
modeling with modes of harmonic degree higher than 3 
necessary. 

Daszynska-Daszkiewicz et al. (2002) developed a new 
method for the determination of I and / from photometric 



amplitudes in different filters and radial velocity ampli- 
tudes. They applied their method on the simultaneously 
obtained photometric and spectroscopic data of FG Vir 
in 2002 (Daszynska-Daszkiewicz et al. 2005). Twelve pul- 
sation modes have been detected spectroscopically in the 
radial velocity curve and all of them have photometric 
counterparts in v and y (Breger et al. 2005). A comparison 
of the empirical and theoretical / values of these modes, 
adopting different descriptions of the mixing length pa- 
rameter a, yielded the best agreement for < a < 0.5, 
meaning that convection is relatively inefficient in FG Vir. 

3. Observations and data reduction 

3.1. New measurements 

During 2002 February to June we carried out an observing 
campaign for the 5 Scuti star FG Vir in the framework 
of the Delta Scuti Network (Rodler 2003, Zima 2002). 
Photometric as well as spectroscopic time-series observa- 
tions were acquired in order to enable the detection and 
identification of more pulsation frequencies, pinpoint stel- 
lar parameters like the inclination and the equatorial rota- 
tional velocity, and consequently improve stellar pulsation 
models. The combination of simultaneous photometry and 
spectroscopy enabled mode identification methods, which 
make use of flux and velocity variations and thus pro- 
vide a test for the efficiency of convection (Daszynska- 
Daszkiewicz et al. 2005). 

To avoid aliasing effects caused by observing gaps dur- 
ing daylight and unpassable weather, multi-site observa- 
tions from different continents were necessary Daily alias- 
ing, which creates spurious peaks in a Fourier spectrum 
spaced by 1 d" 1 , can seriously prevent the detection of 
intrinsic signal. Since FG Vir is located near the celestial 
equator (S = —5°) observations from both hemispheres 
were possible. We selected telescopes of the 2-meter class 
at sites in Germany, France, USA, Chile, South Africa and 
New Zealand which permitted data acquisition of high- 
dispersion (R > 30000) spectra and of high signal-to- noise 
(S/N > 100) covering a large wavelength range. The cho- 
sen integration times are between 7 and 10 minutes, which 
is about 20 % of the shortest detected period of FG Vir 
and 6 % of the dominant frequency at 12.716 d _1 . The ob- 
servations at ESO were made with a shorter integration 
time of 5 minutes to explore the higher frequency domain. 

Spectroscopy of FG Vir was carried out at the follow- 
ing six sites: 

1. Measurements at the Mt. John University Observatory 
(New Zealand) were carried out with the High 
Efficiency and Resolution Canterbury University Large 
Echelle Spectrograph (HERCULES) attached to the 
lm telescope. The wavelength range was 4585- 
5837 A and the resolution R=35000. The integra- 
tion time was 7 minutes. The observers were Duncan 
Wright, John Bentley and Peter Cottrcll. 

2. The South African Astronomical Observatory (SAAO, 
South Africa) data were acquired with the GIRAFFE 
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spectrograph attached to the 1.9m telescope. The 
wavelength range was 4400 to 7000 A at a resolution 
of ~ 32000. The integration time was 7 minutes. The 
observer was Wolfgang Zima. 

3. The observations at the 2.1m telescope of the 
McDonald Observatory (Texas, USA) were made with 
the Sandiford Cassegrain Echelle Spectrograph at a 
resolution of R=60000 in a wavelength range of 4595- 
5194 A. The integration time was 7 minutes. The ob- 
server was Ulrike Heiter. 

4. The measurements at Tautenburg Observatorium 
(Germany) were obtained at the 2m telescope 
equipped with a Coudc Echelle spectrograph with a 
resolution of R=63000 and in a wavelength range of 
4708-7362 A. The integration time was 7 minutes. The 
observer was Holger Lehmann. 

5. At the European Southern Observatory (ESO, Chile) 
the 2.2m telescope was used with the Fiber- fed 
Extended Range Optical Spectrograph (FEROS). The 
wavelength range was 4379-6549 A at a resolution of 
R=48000 and an integration time of 5 minutes. The 
observer was Ennio Poretti. 

6. Observations at the Observatoire de Haute-Provence 
(OHP, France) were obtained using the AURELIE 
spectrograph attached to the 1.52m telescope. The 
spectral domain covers the range 4470-4540 A with 
a resolution power R=55000. The integration time 
was typically 10 minutes. The observer was Philippe 
Mathias. 



A journal of the observations is given in Tabled 



3.2. Data reduction 

The data from SAAO, ESO, McDonald and OHP were 
reduced with standard IRAF (Image Reduction and 
Analysis Facility) routines. We applied bias, flatfield and 
stray-light corrections to the spectra. Every single order 
was extracted and a wavelength calibration by utiliza- 
tion of the ThAr comparison spectra was applied. The 
reduction of the Tautenburg data was carried out in the 
same way as described above with an additional radial ve- 
locity zero point correction using telluric lines. The data 
from Mt. John Observatory were reduced with a standard 
reduction package (Skuljan et al. 2000). An adjustment 
of the wavelength to the heliocentric frame was carried 
out and the times of observation were transferred to the 
Heliocentric Julian Date. Finally, we normalized the spec- 
tra to the pseudo continuum of each order. 

An example of the observed line profile variability 
(LPV) of FG Vir is given in Fig. ^ It shows the inten- 
sity deviation from the mean for a set of high S /N spectra 
in the range 5360 to 5375 A taken during one observation 
night at ESO. The synchronous movement of the absorp- 
tion lines is mainly due to the dominant mode at 12.7 d . 



Table 1. List of spectroscopic observations during the 
2002 FG Vir campaign. 
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Fig. 1. Profile variations of the absorption lines between 
5362 and 5373 A during one night of observation at ESO. 
The residuals after subtraction of the mean profile are 
shown in gray scales with the scale on the right-hand side 
(continuum units). 
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3.3. Selection of lines 

The most important criteria for the selection of the ab- 
sorption lines to be analyzed were identical temporal be- 
havior and the avoidance of blending. 

We calculated model atmospheres and synthetic spec- 
tra using stellar parameters and abundances derived by 
Mittermayer & Weiss (2003), hereafter MW03, to sup- 
port the search for unblended lines. Unfortunately, we did 
not find a single line which met the requirements in the 
overlapping wavelength region of all six observatories. The 
range below 4700 A was omitted for the analysis, due to 
the rapidly decreasing S/N of the SAAO data towards the 
blue. Thus, we decided to omit the McDonald and OHP 
data for the overall frequency analysis and mode identifi- 
cation in favor of more suitable lines. This resulted in a 
loss of potential frequency resolution due to the shorter 
time base, but enabled us to study line profiles of higher 
quality. A separate analysis of the LPV of the OHP spec- 
tra is reported at the end of Section Q] 

We selected four unblended absorption lines shown 
in Fig. for the mode identification. In the lower panel 
the mean spectrum (solid line) is plotted together with a 
synthetic spectrum (dotted line) calculated from a static 
model atmosphere 1 adopting the following input param- 
eters: T e g —7425 K, log g=3.9, micro turbulence velocity 
w m icro=3.9 kms -1 and ?;sinz=21.3 kms -1 (from MW03). 
The upper panel shows the standard deviation from the 
mean. Due to the pulsationally induced intensity varia- 
tions, the deviation is larger at the position of the lines 
than at the continuum. 

The following lines were chosen: Fei A5324.179 A, 
Fen A5362.869 A, Fei A5367.467 A and Fei A5383.369 A. 
The wavelength values of the lines were taken from the 
Vienna Atomic Line Database (VALD) and chosen as zero 
points in the conversion of Angstrom to kms -1 . The data 
were weighted according to their S/N measured in two 
continuum-regions between 5310 and 5400 A. The mea- 
surements taken at ESO, Mt. John and Tautenburg all 
have a similar average S/N of about 200, whereas the qual- 
ity of the SAAO observations is much lower with about 
80. 

We examined the temporal behavior of the four lines 
separately by means of a period analysis of their equivalent 
width variations. An identical temporal behavior is espe- 
cially important when combining line profiles to a single 
profile since deviations would distort the combined profile. 
No significant differences in the radial velocity amplitudes 
and equivalent width variations of the four lines could be 
detected. Therefore, it was safe to combine them by com- 
puting a mean profile after the wavelength values of every 
line were converted to the Doppler scale centered on the 
values given by VALD. The mean S/N of this combined 

1 The synthetic spectrum was computed with the program 
SYNTH (Piskunov 1992) using parameters of the spectral 
lines from VALD (Piskunov et al. 1995, Kupka et al. 1999, 
Ryabchikova et al. 1999). The model atmosphere was calcu- 
lated using the ATLAS9 code (Kurucz 1993). 




i i i i i i i i 

5320 5340 5360 5380 

wavelength [Angstrom] 

Fig. 2. Spectrum between 5320 and 5400 A. The lower 
panel shows the mean spectrum (solid line) of the four 
selected observatories and lines, which were chosen for a 
closer analysis are marked. A synthetic spectrum com- 
puted from a static model atmosphere (MW03) is shown 
as dotted line. The standard deviation from the mean is 
plotted in the top panel. 

profile is 240, which is an appropriate value for the appli- 
cation of the FPF method. 

Finally, the resulting combined profiles were deconvo- 
lutcd with the instrumental profiles, derived for the data 
set of each observatory separately. Such an approach was 
necessary in order to eliminate different line broadenings 
due to the different instrumental profiles. 

4. Period analysis 

We examined the pulsational content of the combined line 
profile by means of a separate analysis of the radial veloc- 
ity and pixel- by-pixel variations. 

Our approach for the detection of the frequencies was 
carried out following the recursive procedure described 
below. At each step a Fourier spectrum of the original 
or prewhitened data was calculated. The peaks with the 
highest amplitudes were selected and examined for their 
statistical significance. From all detected frequencies, a 
multi-periodic least-squares fit was calculated and the 
data were prewhitened with the derived fit. This cycle was 
performed until the amplitude spectrum was void of statis- 
tically significant periodicities. For our data, we adopted 
the criterion that a peak is statistically significant if its 
S/N is above 4.0 in amplitude (Breger et al. 1993). 

The left hand side of Fig. [3] shows the Fourier analy- 
sis of the radial velocity variations of FG Vir. According 
to the spectral window in the top panel, the daily alias- 
ing is low, giving the 1 d _1 -aliasing peaks only 50 % the 
height of the central peak. The raw Fourier spectrum is 
dominated by the frequency fi — 12.716 d -1 also showing 
up in the photometric data as the mode of highest ampli- 
tude. The order of the modes' radial velocity amplitudes 
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Fig. 3. Fourier analysis of the radial velocity variations (left panels) and. the mean Fourier spectrum across the line 
(right panel) of the combined line profile. From top to bottom, every panel shows an amplitude spectrum prewhitened 
by all frequencies marked in the panels above. The significance limit in the bottom panel is calculated according to 
4*noise over a 10 d _1 -range. Note that the scale of the ordinate is changed in the different panels. 



is very similar to the photometric amplitudes. This is not 
really a surprise, since both values are integrated across 
the stellar disk and the detection of low-degree modes is 
favored. During each iteration of the frequency search we 
replaced the frequency values of the detected modes with 
the more accurately known values from the 2002-2004 pho- 
tometry (Breger et al. 2005). Consequently, during the 
least-squares fit the frequency values were fixed and only 
the zero point, amplitude and phase were improved. 

In TableElwe list all 15 statistically significant periodic 
terms between 9 and 34 d _1 detected in this approach. 
Note that in the following we assign to spectroscopically 
detected terms the prefix v, whereas the modes known 
from photometry are denoted with the prefix / as indi- 
cated by Breger et al. (2005). We report the radial velocity 
observations together with a 12-frequency fit in Fig.f^ 

Eleven detected modes correspond to the strongest 
photometric components having y-amplitudes between 22 
and 1 mmag. There are two modes, /12 = 23.397 d -1 and 



fio — 24.194 d -1 having photometric y-amplitudes above 
1 mmag, which do not appear in our frequency so- 
lution. Both modes are components of close frequency 
pairs of which one component has been detected spectro- 
scopically. The frequency resolution of the spectroscopic 
data of 0.019 d _1 does not permit a separation of the 
first pair, but the separation of the two close modes at 
24.2 d _1 should be sufficient. Indeed, the highest remain- 
ing peak of the residuals after subtraction of the 12 modes 
is at 25.21 d _1 , which turns out to be an alias of /io- The 
spectroscopic detection of — 20.83 d _1 , which has a 
photometric y-amplitude of only 0.33 mmag is remarkable. 
No harmonics or combination frequencies could be found. 

We compared the highest remaining peaks, which are 
all below the significance limit, after prewhitening of all 
detected modes with the frequencies detected by Breger 
et al. (2005). The residual power spectra are displayed in 
Fig. The only components appearing also in the 79 pho- 
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Fig. 4. Observations and 12f-fit of the radial velocity variations of the averaged line profile after shift to the heliocentric 
frame. (+ Mt. John, x SAAO, • Tautenburg, A ESO) 



tometric frequencies are /io and = 34.39 d _1 at very 
low S/N. Other residual peaks are apparently not related 
to detected photometric frequencies. There remains an in- 
creased noise level in some frequency domains, especially 
between 20 and 30 d _1 where photometry has revealed 
many additional frequencies. 

The frequency analysis of the pixel-by-pixel variations 
shows a slightly different picture. The frequency spectra 
at subsequent stages of prewhitening are reported in the 
right hand side of Fig- El Here, the mean of every Fourier 
spectrum computed for each wavelength bin across the 
profile is displayed. Since here, in contrast to the Fourier 
analysis of the radial velocity, the phase information across 
the line is neglected, high degree modes can be detected. 

The frequency spectrum is again dominated by vi . The 
subsequent modes have relatively low radial velocity am- 
plitudes. The term ^13 = 12.79 d _1 , which was not de- 
tected in the radial velocity variations, shows up here as 
the fourth highest peak. Intriguing is the detection of the 
combination mode ^15 = 33.55 d -1 = v\ + vyi- This is 
the only term not found among the 79 photometrically 
detected frequencies and it is the only spectroscopically 
revealed combination frequency. After prewhitening the 
data for 14 frequencies, none of the remaining peaks can 
be identified with a photometric counterpart. 

We determined the mean radial velocity of FG Vir rel- 
ative to the heliocentric frame by cross correlating the 
mean of the observed spectrum with a synthetic spectrum 
between 5310 and 5390 A. Using this method the velocity 
information of all lines in that range can be utilized. The 
location of the main peak in the cross correlation function 



Table 2. Frequencies detected by Fourier analysis of the 
radial velocity (RV) and pixel-by-pixel (PbP) variations. 
The frequency values are taken from photometry (Breger 
et al. 2005). The amplitudes are calculated from a multi- 
periodic least-squares fit to the data. The PbP amplitude 
is the integral of the amplitude across the line. The photo- 
metric amplitudes are derived from a 60 frequency least- 
squares fit to the 2002 data. 



id 


freq. 


RV amp. 


PbP amp. 


y amp. 


v amp. 




d- 1 


kms^ 1 


kms^ 1 


mmag 


mmag 






±0.03 


±0.02 


±0.07 


±0.08 




9.199 


0.28 


0.33 


2.73 


3.80 




9.656 


0.28 


0.20 


3.62 


5.14 


vz 


12.154 


0.36 


0.21 


4.10 


6.08 


V\ 


12.716 


2.06 


1.16 


22.10 


32.13 




12.794 




0.29 


0.68 


0.82 


V\\ 


16.071 


0.16 




1.23 


1.61 


"9 


19.227 


0.18 


0.29 


1.73 


2.52 


V7 


19.867 


0.31 


0.20 


1.96 


2.79 


V10 


20.287 


0.23 


0.35 


1.18 


1.84 


V±2 


20.834 


0.15 


0.27 


0.33 


0.49 


VA 


21.051 


0.38 


0.23 


3.07 


4.32 


Vb 


23.403 


0.36 


0.25 


3.98 


5.33 


V-2 


24.227 


0.43 


0.25 


4.26 


5.65 


via 


25.433 




0.16 


0.98 


1.34 




33.551 




0.18 







provides the shift between these two spectra. Thereby we 
obtained a value of ±14.4±0.2 kms -1 . 

The measurements at ESO have a short integration 
time of 5 minutes reducing the effects of phase smear- 
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Fig. 5. Residual power of the radial velocity between 5 
and 50 d _1 . Left panel: The residuals after prewhitening 
12 frequencies. Two of the highest remaining peaks are 
identified with modes known from photometry (the term 
at 25.2 d _1 is a one-day alias of /io). Right panel: After 
additional prewhitening of these two frequencies, none of 
the remaining peaks can be identified with a photometric 
counterpart. Especially in the regions around 14 d _1 , from 
20 to 28 d _1 , around 35 d _1 and 43 d _1 the noise level is 
noticeably increased. A similar situation was observed in 
the photoelectric data from 1995 (Breger et al. 1998). 



ing for higher frequencies. We carried out a separate fre- 
quency analysis of the radial velocities and pixel-by-pixel 
variations. No additional statistically significant frequen- 
cies above ^15 = 33.55 d _1 could be detected spectroscopi- 
cally from these data. More multi-site observations having 
such short integration time would be required to permit an 
improved frequency analysis of this high frequency range. 

Comparing our frequency analysis to Mantegazza & 
Porctti (2002), hereafter MP02, the striking difference to 
their results is the radial velocity amplitude of 1/1. A re- 
analysis of the 1995 data 2 has revealed some inconsisten- 
cies in the determination of the radial velocity values as 
obtained from the cross correlation of the normalized spec- 
tra with the average spectrum. Indeed, when calculating 
the radial velocity amplitudes as first moment (i.e., the 
line barycentre) the agreement with the values listed in 



Table 2 is excellent 3 . Therefore, there is no reason to sus- 
pect the variability of the radial velocity amplitude of the 
v\ term, also considering that its photometric amplitude 
is stable. Moreover, we note that the analysis of the 1992 
data (Mantegazza et al. 1994) supplied very similar radial 
velocity amplitudes. 

Interestingly, MP02 were able to detect three photo- 
metric frequencies at 34.12, 21.23 and 21.55 d _1 , which do 
not show up in our data. An explanation might be that 
they used deeper absorption lines in their analysis result- 
ing in a higher ratio of pixel-by-pixel amplitudes to the 
observational noise. 

A separate frequency analysis of the radial velocities 
and line profile variations was also undertaken for the 
OHP spectra using the mean of the two metallic lines Ti 11 
4501 and Fe 11 4508 A. Radial velocities were computed us- 
ing two methods: a Gaussian fit, and the first normalized 
moment. 

A Fourier analysis clearly shows the v\ frequency, 
which has an amplitude of 2.1 kms -1 . After prewhiten- 
ing with ui, many peaks appear, the dominant one being 
around 21.2 /d (alias of 1/2). Many significant peaks are 
detected: 12.15 (fa), 9.66 (/s), but they are polluted by 
aliasing problems due to the mono-site situation. Table 
3 lists the radial velocity amplitudes of the OHP spectra 
computed by a simultaneous multi-periodic least-squares 
fit taking the frequencies from Table 2. 

Table 3. Radial velocity amplitudes of the OHP spectra 
computed from a multi-periodic least-squares fit using the 
frequencies listed in Table 2. 



frequency 


RV amplitude 


d- 1 


kms -1 


12.716 


2.10 


24.227 


0.33 


12.154 


0.51 


21.051 


0.52 


23.403 


0.31 


9.199 


0.18 


19.867 


0.13 


9.656 


0.25 


19.227 


0.37 


20.287 


0.23 


16.071 


0.22 


20.834 


0.21 


12.794 


0.19 


25.433 


0.19 


33.551 


0.04 



In 1999, single-site high-resolution Echellc spec- 
troscopy of FG Vir was carried out during the nights 
of May 3, 6, and 7 at McDonald Observatory by T. 



Kindly provided by L. Mantegazza 



3 The new radial velocity amplitudes for the 10 terms so- 
lution provided by MP02 in their Table 1 are, from top to 
bottom: 2.10, 0.51, 0.31, 0.53, 0.34, 0.41, 0.35, 0.34, 0.38 and 
0.13 kms" 1 
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J. Montemayor. Fburtytwo spectra between 5500 and 
6000 A, having an average S/N of 220, were acquired. We 
used these data to check the stability of the radial velocity 
amplitude of v\ and the mean radial velocity relative to 
the heliocentric frame. 

Only v\ could be detected in the variations of the radial 
velocity and pixel intensities of the Fe I 5525 A line. Other 
periodicities could not be resolved due to the small number 
of data points. We derived a radial velocity amplitude of 
1.81±0.2 kms -1 , which is consistent with the value derived 
in 2002. 

The mean radial velocity, + 14.7±0.5 kms^ 1 , was de- 
rived from lines in the range between 5520 and 6002 A by 
the same procedure as described above for the 2002 data. 
Therefore, we can say that no significant radial velocity 
shift occurred between the three years. 

Also considering that a new analysis of the 1992 and 
the 1995 data supplied mean radial velocity values of 
14.4 ± 0.3 and 15.4 ± 0.2 kms -1 , respectively, we can say 
that no significant radial velocity shift occured between 
1992 and 2002. 



of the MM and by Schrijvers et al. (1997). The intrinsic 
unbroadened line profile was assumed to be a Gaussian. 
The broadened profile was computed by summing over a 
7200-element grid approximating the visible distorted stel- 
lar surface. A quadratic limb darkening law was adopted. 
For more details on the computation of the theoretical 
LPV refer to Appendix A of Paper I. 

The observed uncertainties cr(A\) and cr(4>\) were com- 
puted from the standard deviation of the pixels across the 
line after prcwhitcning of all significant frequencies. The 
behavior of the observed zero point, amplitude, and phase 
of all 15 detected terms is shown in Figure^! The observed 
values are denoted by the crosses with error bars, whereas 
the lines denote fits described below. The top-most panel 
is the zero point profile, which is identical for all frequen- 
cies. For each frequency's diagram, the upper panel shows 
the amplitude in units relative to the continuum and the 
panel below shows the phase distribution expressed in ra- 
dians. 

A visual inspection on these diagrams already reveals 
a lot of interesting information about the mode properties. 



4.1. Mode identification by means of the FPF method 

For pinpointing the properties of the pulsation modes we 
applied the FPF method. By means of this method, the 
observed amplitude and phase values across the profile 
are fitted with theoretical values derived from synthetic 
line profiles. The fitting is carried out by applying genetic 
optimization routines in a large parameter space. The re- 
duced xt is calculated from complex amplitudes in order 
to combine amplitude and phase information in the fol- 
lowing way 



xl = 



1 



2n x 



1=1 



'R.i 



(1) 



where n\ is the number of pixels across the profile, m 
is the number of free parameters, A° and A 1 denote ob- 
servationally and theoretically determined values, respec- 
tively, Afj — A\ cos 4>\ and Aj = A\ sin <j>\ are the real and 
imaginary part of the complex amplitude, and a 2 is the 
observed variance. Eq. JIJ can be easily modified if the ob- 
served amplitude and phase from photometric passbands 
should be included for the calculation of xl- 

Since the amplitude and phase of a given wavelength 
bin are treated as independent variables, the variances are 
calculated from 



2 

a R,X 



'I.X 



-- <t(A x ) 2 cos 2 . 
a (A x ) 2 sin 2 4 



) 2 A 2 sin 2 , 



b x ) 2 A 2 cos 2 , 



(2) 
(3) 



We describe the theoretical stellar surface displace- 
ment field as superposition of spherical harmonics assum- 
ing linear pulsation theory. The description is limited to 
a slowly rotating system taking into account toroidal mo- 
tion caused by the Coriolis force. Such a description has 
also been applied by Aerts & Waelkens (1993) in the case 



— The dominant mode i^i = 12.72 d _1 shows clearly 
an axisymmctric behavior (m = 0), namely a double- 
peaked amplitude distribution with zero amplitude at 
the center of the profile accompanied by a phase shift 
of 7r at the same position. The difference of the height 
of the two amplitude maxima is intrinsic to the LPV. 
It is caused by the variable equivalent width due to the 
temperature variations. 

— The double-peaked amplitude distribution of five other 
modes, namely v$ — 9.66 d _1 , u-j — 19.87 d _1 , 



21.05 d~\ v 5 = 23.40 d _1 , and v 2 = 24.23 d _1 , also 
shows the characteristics of axisymmetric modes. Here, 
differences in the heights of the two maxima might 
be due to observed uncertainties considering the error 
bars. The phase-jump of 7r at the minimum amplitude 
is also consistent with modes of m = 0. 

— This high number of axisymmetric modes is in agree- 
ment with the observations reported by MP02, who 
have identified 8 out of 10 detected modes as m = 0. 
The reason for this might be a relatively low incli- 
nation angle combined with the low v sin i of FG Vir 
favoring the detection of axisymmetric modes, which 
are best seen pole-on. Since the rotationally broadened 
line profile represents an one-dimensional projection of 
the stellar surface, higher values of v sin iresult in spa- 
tially better resolution. Therefore, high-degree modes 
are easier detectable at high v sin iand FG Vir at about 
vshii=21 kms _1 is not expected to show spectroscop- 
ically high-degree modes, just as it is observed. 

— At least five modes are non- axisymmetric: v§ — 
9.20 d" 1 , v 9 = 19.23 d~\ v%o = 20.29 d" 1 , u 12 = 
20.83 d" 1 , and i/ 13 = 12.79 d -1 . This is deduced from 
their amplitude distribution, which has a maximum 
close to the line center, combined with a gradual phase 
change across the line of more than it and therefore an 
azimuthal order m > 0. i/q, v$ and v\ 2 appear to be 
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retrograde modes (m > 0) due to the gradual increase 
in phase from blue to red, whereas and v\q are 
apparently prograde modes. 

The diagrams of = 12.15 d" 1 do not permit 
a conclusion about its character. Photometric mode 
identification (Breger et al. 1999) and the results by 
Daszynska-Daszkiewicz et al (2005), hereafter refered 
to as DD05, strongly point towards the identification 
as fundamental radial pulsation mode. From our de- 
rived amplitudes and phases across the line, we can- 
not confirm these results. A reason might be that in 
our spectroscopic data it is impossible to separate 
1/3 from a photometrically detected close frequency 
at 12.162 d -1 . An interaction of these two modes in 
Fourier space could prevent a successful mode identi- 
fication. 

Due to the low S/N of v\\ = 16.07 d _1 we exclude 
this mode from our mode identification. This includes 
also the harmonic and combination terms, v\± and 
^15 respectively, since these are not intrinsic pulsation 
modes. 



E 0.9 - 




+ observational zero-point profile 
— fit(vsini=21.6km/s, EW=8.0 km/s,a=1 1 .2 km/s) 

i l i l i I i I 

) -20 -10 10 

velocity [km/s] 



Fig. 7. Observed zero point profile and best fit by means 
of a synthetic rotationally broadened profile. 



Our modeled line profiles depend on the stellar pa- 
rameters: mass M, T e g, logg, limb darkening, and on the 
flux derivatives, olt and a g . We adopted the mean stellar 
parameters derived by DD05. Photometric calibrations of 
Kurucz models as well as Vienna atmospheres yield con- 
sistent values of log T cfr = 3.87±0.01 and log 5 =4.0±0.1. 
DD05 determined an improved value of T e ff by fitting 
the amplitude and phase of the dominant mode v\ to 
the observations using flux derivatives determined from 
NEMO.2003 models (Nendwich et al. 2004). The best fit 
was achieved for a model of solar metal abundance and a 
micro-turbulence velocity of £ — 4 kms -1 , yielding a value 
of log T c g =3.86 and logg=4.0. Our flux derivatives ax 
and a g were computed from NEMO.2003 models choosing 
the previously determined values of T e ff and logg. 

Our limb-darkening coefficients were taken from tab- 
ular data provided by Barban et al. (2003). We adopted 



the values for the quadratic limb-darkening law for the 
Stromgren filters v and y. We excluded the wings of the 
line from our analysis due to the fact that there the ro- 
tationally broadened profile is dominated by the Lorentz 
profile deviating from the adopted Gaussian intrinsic pro- 
file. Therefore, our fit covers only the range between -26 
and 26 kms -1 . 

The projected rotational velocity wsini, the width of 
the intrinsic profile a, and the equivalent width W were 
computed from a least-squares fit of a rotationally broad- 
ened synthetic profile to the zero point profile (see Fig.EJ). 
We obtained the following values: v sin?=21.6±0.3 kms -1 , 
VF=8.02±0.05 kms" 1 , a = 11.2 ± 0.3 kms" 1 . The zero 
point profile is slightly broadened due to the LPV, mainly 
by the contribution of v\, which has RV/usini=0.1. The 
resulting overestimation of usini is only about 1%. 



4.2. Results of the mode identification by neglecting 
equivalent width variations 

We started our mode identification by neglecting the 
equivalent width variations of the spectral line due to 
temperature variations. This was achieved by setting 
a\y = 0. A genetic optimization (for details see Paper 
I) with the following free parameters [min, max; step] 
was adopted: (—[0, 4; 1], m—[-£, £; 1], the relative dis- 
placement amplitude a s =[0.00005, 0.003; 0.00005], i=[5, 
90; 5]°, usini=[21.3, 21.9; 0.1] kms" 1 , cr=[10.9, 11.5; 
0.1] kms -1 . Since we know the equivalent width of the 
line with much better precision than wsini and <r, we 
fixed it during the optimization to 8.02 kms -1 . The theo- 
retical equatorial break-up velocity of FG Vir is at about 



= 210 kms -1 (Collins & Harrington 1966), which leads 
to a minimum possible inclination angle of about 6°. It 
therefore is reasonable to set the lower inclination limit to 
this value. For each pulsation mode we calculated approx- 
imately 40000 models during the process of optimization. 

According to this identification, six modes are axisym- 
metric, four are retrograde with m = 1 and two modes 
are prograde with m = — 1 and m = — 2 respectively. The 
inclination of the best fitting models covers almost the 
whole tested range from 10 to 85°. 

Only four of the twelve tested modes, ^5, vq, and 
^13 can be fitted satisfactory in the adopted confidence 
limit. The \i of v\ is especially high (xi — 10.72), which 
is mainly due to the fact that this mode clearly shows 
equivalent width variations which are ignored here. We 
will focus on this problem in detail in the next subsection, 
where we discuss the fits when equivalent width variations 
are included. Nevertheless, for nine modes the three best 
solutions yield identical m values. The fact that m is much 
better constrained than I is consistent with the results of 
our experimental tests with synthetic data (see Paper I). 

There are some modes, especially 1^7 = 19.87 d -1 , 
v\q = 20.29 d _1 and 1^2 = 24.23 d _1 , showing a signifi- 
cant velocity shift between the observed and modeled am- 
plitude distribution that might be connected to equiva- 
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Table 4. Mode parameters derived from the application of the FPF method by including equivalent width variations. 
For each pulsation mode, the five best solutions are shown. a s is the relative displacement amplitude; i is the stellar 
inclination angle in degrees; vs'mi is the projected rotational velocity, a is the width of the intrinsic profile, both 
expressed in kms -1 ; |/| and ipf are the ratio and the phase shift, respectively, between flux and radius variations. 
aw describes the dependence of the equivalent width with respect to the temperature. The 95 % significance limit is 
at xi — 1-37. Due to spherical symmetry we have not indicated the inclination angle for radial modes. 
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lent width variations of these terms. Another explanation 
considers photometrically detected close frequency terms 
which distort the Fourier parameters across the line. 

4.3. Results of the mode identification by including 
equivalent width and light variations 

In an improved approach we carried out a mode identi- 
fication by also including the equivalent width variations 
in a parameterized form caused by the pulsational tem- 
perature variations. The complex parameter /, the ratio 
of the stellar flux to radius variations, governs the rela- 
tive temperature change and therefore also the equivalent 
width change of the line. This gives us the opportunity to 
derive / and its phase shift ipf empirically from the ob- 
served equivalent width variations. The parameter space 
for the optimization was extended by varying the follow- 
ing values: |/| [0,22;0.5], ip f [0, 2tt;0.05], the change of the 



equivalent width with respect to the temperature change 
% [-5, 5; 1] (see Paper I for details). 

We report the best fits for each mode in Tables 0] and 
the corresponding amplitude/phase diagrams in Figure |SJ 
The xi value of five modes, 1^4, ^5, vq, 1/9, and 1^13 reaches 
below the 95 %-significance limit of 1.37. Now, addition- 
ally iAj = 19.23 d _1 can be satisfactorily modeled as a 
retrograde mode of m = 1 and I = 1 or 2. For all modes, 
the minimum xi is lower than in the first approach that 
did not include equivalent width variations. 

Most importantly, the fit for the dominant mode V\ 
is much better now and points towards a dipole mode, 
consistent with previous mode identifications from other 
authors (Breger et al. 1999, Mantegazza & Poretti 2002). 
From the fit to this mode an inclination of 35° is derived. 
But, as shown in Paper I, axisymmetric modes are inap- 
propriate for determining the inclination since almost any 
inclination value can be obtained by scaling the ampli- 
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tudc. It is much better to rely on the non-axisymmetric 
modes like vq, and ^13, which all point towards a low 
inclination value between 10 and 25°. We will explore this 
in more detail in the next subsection, where an improved 
value of i is computed by simultaneously fitting these three 
modes. 

Formally, we can only fit five pulsation modes satisfac- 
torily and other modes can not be reproduced within the 
adopted significance limit. But also for the latter modes, 
some I and m must be correct. Apparently, their Fourier 
diagrams are distorted in such a way that we cannot model 
them. 

In most cases, the fact that the modeled Fourier val- 
ues fail to fit the observations is due to a blue or red shift 
of the observed amplitude and phase values. This is best 
seen for v-i and v\^,. Identical shifts have also been seen 
in the Fourier parameters of single absorption lines so the 
combination of the lines can be ruled out as an explana- 
tion. Also the amplitude/phase diagrams given by MP02 
indicate such "velocity shifts" . 

Close unresolved frequency pairs can cause such a sit- 
uation if spectroscopically only one component is visi- 
ble (see Paper I). Looking at the 79 photometrically de- 
tected frequencies by Breger et al. (2005) we see that ^2, 
v 3i v 5> V 7 A and vw indeed have close components, which 
might distort the amplitude and phase distribution across 
the line if they are not taken into account in the multi- 
frequency least-squares solution. 

We calculated new Fourier parameters from a least- 
squares fit containing 23 frequencies in the solution. All 
photometrically detected frequencies that are closer than 
0.03 d _1 to spectroscopically found frequencies were in- 
cluded. The results show that only the amplitude/phase 
distribution of 1/3 and v\q exhibit a significant change. By 
means of the FPF method ^3 then can be identified as a 
radial or dipole mode, as also indicated from photomet- 
ric mode identification (Breger et al. 1998) and by DD05, 
whereas v\q has the best solution at m — —2 in a close 
frequency solution. Unfortunately, the values of \t are 
not improved compared to the 14f-solution. An additional 
problem is that a slight change of the close frequency's 
value has a large impact on the amplitude/phase values 
across the line of the mode under consideration. 

Both l/l and ipf are relatively poorly determined for 
most modes. The xt diagrams indicate a wide range of 
possible / values. This is also true for aw, which is di- 
rectly dependent on /. We have shown that very high 
S/N data are needed to obtain reliable values of / from 
the equivalent width variations alone. 

Taking additional photometric measurements into ac- 
count significantly improves the quality of the derived /. 
DD05 have shown that it is essential that both spectro- 
scopic and photometric data are obtained in the same year 
due to the fact that slight year-to-year changes of ampli- 
tude and phase may exist. These changes might be due to 




4 It could not yet be proven if this mode has a close compo- 
nent or if its amplitude is variable (Breger et al. 2005) 



16 20 
frequency [c/d] 

Fig. 8. Comparison of the empirically derived real (top 
panel) and imaginary (bottom panel) part of the complex 
parameter / with theoretical values for different mixing 
length parameters a. The empirical values are shown as 
black dots with error bars. Each line denotes theoretical 
values for different values of a. 



intrinsic amplitude variability of the star and/or to unde- 
tected low-amplitude frequencies. Therefore, we adopted 
all 79 detected frequencies to a least-squares fit by only 
taking into account the photometric data from 2002. We 
adopted the FPF method by simultaneously fitting the 
Fourier parameters of each mode and its photometric am- 
plitude/phase in v and y. Compared to the previous fits, 
there is no significant difference in the rank of the best 
identifications of £ and m for all modes, however, for most 
of the modes |/|, %pf and aw are better constrained now. 

It has been shown by Daszynska-Daszkiewicz et al. 
(2003) that the computed values of / are sensitive to the 
description of convective transport. In 5 Set stars like FG 
Vir two unconnected convective zones associated with the 
H and Hell ionization zones exist. The energy transport 
in these convective layers governs the flux variations at 
the photosphere. A commonly accepted approach to de- 
scribe the complex convective process is given by the mix- 
ing length theory. 

In Figure|Hlwe present our empirically derived real and 
imaginary part of / together with theoretical values for 
different mixing length parameters a. Note that the sign 
of our imaginary // is inverted compared to DD05 due to 
the fact that we adopt different conventions of the time 
dependence e tuJt of the spherical harmonics. The theoret- 
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ical values were computed with OPAL opacities (Iglesias 
& Rogers 1996) under the assumption of standard Pop I 
composition (Y = 0.28, Z = 0.02). DD05 have shown that 
the choice of Z has little effect on the computed / values, 
whereas an increase of Y to 0.33 causes a shift of f R and // 
in the order of 2. The observations agree well with models 
having a low mixing length coefficient a < 0.5 meaning 
that convection is relatively ineffective in FG Vir. This is 
in good agreement with independently obtained results by 
DD05. 

The results of our mode identification by means of 
the FPF method are listed in Table [5j For each mode we 
have indicated the most probable values of £ and to from 
our analysis, together with the independently derived t- 
values from DD05 and m-values from MP02 for compari- 
son. This table emphasizes the need for independently ap- 
plied identification techniques. Whereas I is much better 
constrained from the combined photometric/spectroscopic 
identification technique of DD05, it can not determine 
the value of to. Except for the m-values of the results 
of MP02 are coincident with our identification of the az- 
imuthal order. 

Table 5. Best identification of I and to for FG Vir. 



name 


frequency 


this work 


DD05 


MP02 




d" 1 


I 


m 


I 


m 




9.199 


1, 2, 3 


1 


2 


1 


^8 


9.656 


0, 1, 2 





2 





vz 


12.154 


0, 1, 2 


0, 1 








V\ 


12.716 


0, 1 





1 





V\Z 


12.794 


2, 3, 4 


-2 


1, 2 






19.227 


1, 2 


1 


0, 1, 2 





U 7 


19.867 


0, 1, 2 





2, 1 





V10 


20.287 


1, 2, 3 


-1 


0, 1 




V\2 


20.834 


2, 3, 4 


1 






Vi 


21.051 


0, 1, 2 





1, o 





v 5 


23.403 


2 





2, 1 







24.227 


0, 1 





1 






4.4. An improved determination of the inclination 

For the three non-axisymmetric modes vq, vq and ^13 we 
calculated a simultaneous fit with common values of i, 
usini, o and W. Such an approach is also carried out 
in the case of the moment method by Briquet & Aerts 
(2003). The advantage lies in the simple fact that these 
four parameters are the same for all these modes. A ma- 
jor drawback of such an approach is the large number of 
free parameters, which significantly increases the compu- 
tational time for the determination of the optimum value 
of the inclination. 

We selected the three modes because they can all be 
fitted within the 95 %-significance limit. Moreover, they 
are non-axisymmetric, implying that there is a stronger 
dependence of xt with respect to i than for axisymmetric 
modes (see Paper I). 
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Fig. 9. Determination of the inclination of FG Vir from a 
simultaneous FPF to the three non-axisymmetric modes 
uq, vg and 1/13. Every single symbol represents a computed 
model. The 95 %-confidence limit of xl is indicated as 
solid horizontal line. 



We applied the multi-mode fit in two steps: with and 
without taking into account equivalent width variations. 
We limited the values of / and m to the three best solu- 
tions from the mono-mode fits in order to speed up the 
optimization. 

Figure El shows a clear minimum of the xt values at 
i = 19° ± 5°. The fit is better for the models that include 
the equivalent width variations but the derived inclination 
is the same. Such a small aspect angle was also suspected 
by MP02, who derived i = 15 ± 5 from a line profile fit to 
the variations of five modes. 

Taking into account vsini=21.6±0.3 kms -1 , we derive 
an equatorial rotation velocity of v ro t = 66 ± 16 kms -1 . 
With an assumed radius of 1.58 * 10 6 km (Breger et al. 
1999) this results in a rotational period of fl — 1.74±0.43 d 
and a rotational frequency of 0.57 ± 0.14 d _1 . 

We also see one triplet of rotationally split pulsation 
modes consisting of i>g, v-j and v\§. From such a triplet 
it is possible to derive an approximate rotation frequency 
by considering half of the frequency difference between the 
to = —1 and the m = +1 mode. The computed value of 
0.53 d _1 is very close to the value independently derived 
above. 

5. Conclusions 

We carried out a mode identification for the 5 Set star 
FG Vir using high-resolution time-series spectra. 15 fre- 
quencies could be detected from the radial velocity and 
pixel-by-pixel intensity variations of a profile made from 
the combination of four iron lines. 

We have presented the first successful application of 
the FPF method to line profiles of a 5 Set star. The pul- 
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sational geometry of 12 independent frequencies was an- 
alyzed by means of this method. Seven frequencies were 
identified with axisymmctric modes having I values be- 
tween and 2. Five non-axisymmctric modes with m val- 
ues between -2 and 1 were detected. These modes could 
be modeled with t values as large as 4. The observed 
radial velocity amplitudes range between about 2 and 
0.15 kms" 1 . The best models fit the profiles with 0.17 % 
relative radius variation for the dominant mode v\ and 
down to a value as small as 0.02 %. No high degree modes 
were found mainly due to the low rotational broadening 
of the line. 

By simultaneously fitting amplitudes and phases of ve- 
locity and light variations we computed /-values for eleven 
frequencies. A comparison of these empirical /-values with 
theoretical values computed for different mixing length 
values showed that convection seems to be relatively inef- 
fective in FG Vir. 

The good identifications of the azimuthal order, com- 
bined with the determination of the inclination angle 
i w 19°, will help to improve the treatment of rotational 
effects by seismic models. 

The present fit of the theoretical models to the observa- 
tions cannot be improved significantly by carrying out ad- 
ditional photometric measurements of FG Vir. Hundreds 
of unstable modes with I < 5 are predicted by theory for 
the frequency range between 7 and 40 d _1 (Pamyatnykh, 
priv. comm.). Reducing the observational detection limit 
from 0.2 to 0.1 mmag would yield many more frequen- 
cies but without a knowledge of their pulsational quan- 
tum numbers the number of free degrees is increased with 
every detected frequency. A more promising approach is 
to start the search for a good stellar pulsation model with 
the anchor points given by our mode identification. 
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Fig. 6. Fit of the Fourier parameters across the line profile of the detected pulsation modes taking equivalent width 
variations into account. The top panels show the zero-point profile. For every single frequency, the observed amplitude 
in units of the continuum and the phase distribution in radians is shown together with the two best fitting models. 
For the three frequencies shown in the bottom panels, no fits have been computed due to the low S/N (in the case of 
vn = 16.07 d _1 ) and due to their nature as harmonic (^14 = 25.43 d _1 ) or combination frequency (1/15 = 33.55 d _1 ), 
respectively. 



